Power laws in surface state LDOS oscillations near a step edge 
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In this paper we indicate a general method to calculate the power law that governs how electronic 
LDOS oscillations decay far away from a surface step edge (or any local linear barrier), in the energy 
range when only 2D surface states are relevant. We identify the critical aspects of the 2D surface 
state band structure that contribute to these decaying oscillations and illustrate our derived formula 
with actual examples. 
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I. INTRODUCTION 

For over a couple of decades now Scanning- Tunneling 
Microscopy (STM) experiments have been used to ob- 
serve the effects of perturbations to electronic surface 
states, in the form of atomic defects, corrals and step 
edge^S. The quantum electronic response to defects 
can give us basic information about the band structure 
of the scattered quasiparticles^; it also encodes informa- 
tion about their nature and this is useful when probing 
correlated phase^I^M^. In this paper, we shall consider 
the case of a step edge on a 2D surface and calculate the 
spatial decay of standing waves created in the surface 
LDOS, far away from the step edge. We shall show, by 
a simple process of power counting, that the geometry of 
the constant energy cut of the quasiparticle band struc- 
ture and a qualitative knowledge of the character of the 
quasiparticle wavefunctions provide enough information 
to pin down the power with which the LDOS oscillations 
decay far away from the step edge. 



II. THEORY 

We begin by considering a general band structure for 
the 2D surface states, whose cross-section at the energy 
of observation £^obs is shown in Figure [l] Also shown in 
the figure is the orientation of the surface step edge — 
parallel to the ^/-direction. Because of the conservation of 
momentum parallel to the edge during a scattering pro- 
cess, the incoming and outgoing states must be connected 
by straight lines perpendicular to the step edge. Some 
such processes are also marked in Figure [l] The arrows 
joining the initial and final states denote the wave- vector 
of LDOS oscillations that particular scattering process 
would give rise to. Obtaining the total LDOS involves 
summing up these oscillations. The most coherent con- 
tributions to this sum come from regions where the scat- 
tering wave-vectors change the slowest as we move par- 
allel to the step edge, i.e, changing only the ky of the 
scattering states. We denote the 'identifying' scattering 
wave vector in each such region as the 'characteristic' 
wave vector of that region. 



As a very common example, for a circular constant 
energy cut as in a 2DEG (Figure |2|, the most coherent 
contributions come from the scatterings around the di- 
ameter — the characteristic wave vector in this case is 
thus the diameter, Aq. 



The new electronic LDOS far away from the step edge 
is now provided by (below, 'new' refers to the new energy 
eigenstates while 'init' and 'fin' refer to the initial and 




FIG. 1: (Color online) A constant energy cut of a generic 
2D electronic surface state (SS) band structure, taken at the 
energy Eohs of an STM probe; on the right is shown the ori- 
entation of the surface step edge (or any linear barrier) in 
question. Quasiparticles are scattered 'horizontally', preserv- 
ing ky. The regions where the scattering wave- vectors vary 
the slowest (locally) with ky are shaded - blue for incident 
and pink for reflected states. The 'characteristic' scattering 
wave vectors are also indicated by arrows and numbers. 
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Writing the energy- momentum eigenstates as ipk{x) = 
^^ikx^ where x denotes an 'internal' part involving the 
spin and other internal components, the x-dependent 
part of the LDOS can be summarized as 



M«.^)«2:/*.<'.('.)He[.-fe)(.rx.) 



^iAkxX 



(2) 

The sum is over the various regions of coherent scatter- 
ing, each one corresponding to a characteristic scattering 
vector. The outer limits of these integrals are not impor- 
tant as the oscillations there de-cohere rapidly. po{ky) 
is a DOS factor (it multiplicatively converts the measure 
dky to a product of the length of the band curve enclosed 
between ky and ky + dky and the DOS in that region). 

The x-dependence of a characteristic oscillation far 
away from the step edge may be found from the 
above expression by writing down the lowest order ky- 
dependencies of the relevant quantities near each charac- 
teristic wave-vector {Sky is the /c^y -displacement from the 
associated characteristic wave vector): 
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(3) 



Changing the integration variable 6ky to the variable /i = 
Sk^ X in ([2]) and using ([s]), we obtain our central result 



This exercise provides us with the combination of possi- 
ble power laws we can try to fit actual experimental data 
to, if an idea of the band structure exists. 



III. EXAMPLES 
A. 2DEC^ 

The scattering wave-vector varies slowest around the 
equator (see Figure [2]). Thus, Aq = diameter of circle 
= 2/c^;, Of = /3 = 7 = (assuming, quite reasonably, that 
the reflection amplitude is nonzero and smooth across 
normal scattering). Also, from the geometry of the band, 
we get T] = 2. Using Q, this tells us that: 

Spi.,E)r.'^^^i^^ (Ao = 2M (5) 
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This asymptotic behavior is correct for x ^ (A/c)~^, 
where Ak is the characteristic size of the region in mo- 
mentum space where the scaling laws ([3| hold. The 
power of decay of the oscillations coming from each char- 
acteristic region is thus given by (a + /3 + 7 + 1)/^, which 
may be evaluated using our knowledge of the band struc- 
ture in that region. Note, however, that we haven't been 
able to evaluate the total strength of the scattering pro- 
cess which requires a much more detailed calculation in- 
cluding evaluating the reflection amplitudes themselves. 



FIG. 2: (Color online) The characteristic wave- vector Ao in 
the case of a circular band (for 2DEGs with rotational invari- 
ance) 



B. Strong Topological Insulator (circular band cut) 

1. Generic barrier 

This case is illustrated in Figure [3] and is realized for 
he gapless surface states in Strong Topological Insula- 
tors like Bi2Se3 and Bi2Te3 (at energies near the Dirac 
point). The scattering wave- vector varies slowest around 
the equator (as in the 2DEG case), where Ao=diameter 
of circle, a = 0, /3 = 1 since the reflection coefficient 
changes sign^ as one crosses the diameter/case of normal 
reflection (can be any odd power; should be linear gener- 
ically)^ 7 = 1 because the spins are exactly antiparallel 
for scattering states at the diameter and thus the lowest 
order overlap is linear in 5ky^ and ry = 2 as in the 2DEG 
case. This gives rise to: 
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(Ao = 2kE) (6) 



This result agrees with numerical calculations for partic- 
ular cases of the model describing the step edge^. 




Same contribution 
due toTR 



circle, Qf = 0,/3 = 1,7 = 1 and rj = 2 exactly as argued 
before for the circular STI band. However, the extent of 
this region is very small and the scattering is found to be 
dominated by processes connecting the hexagonal 'cor- 
ners' (marked by bold arrow in Figure |4|, with a charac- 
teristic scattering vector knest^- For the latter case, since 
the spin states have a finite overlap with each other at 
the hexagon corners^^, we have 7 = 0. Also, assuming 
that the refiection coefficient is smooth for the relevant 
scattering processes, /3 = 0. Finally, a = (DOS is finite 
and smooth) and the overwhelmingly 'linear' nature of 
the bands yield 77 = 1. Putting these together, we obtain 
the observed variation^ 
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FIG. 3: (Color online) Constant energy cut of a circular sur- 
face band on a STI surface. The spins, indicated as block 
arrows, are antiparallel for normal scattering (the character- 
istic scattering process for the STI surface state band that is 
circular) - the spin overlap magnitude is thus generically a 
linear function of the angle of incidence. The same may be 
said for the reflection amplitude magnitude which is a certain 
gauge is antisymmetric in the angle of incidence. 



2. 'Perfect' reflection 

The scattering wave-vector varies slowest around the 
equator (as in the 2DEG case), where Ao=diameter of 
circle, a = /3 = 0^ (since the refiection amplitude is con- 
stant in magnitude near normal incidence), 7 = 1 as 
argued for the previous case, and r] = 2. This gives rise 
to: 
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FIG. 4: (Color online) Scattering processes from a step edge 
on Bi2Te3 oriented perpendicular to the FM direction, in the 
energy range where the band exhibits hexagonal warping. The 
scattering vector varies linearly near /cnest , as indicated by the 
angle made by the dotted lines, leading to 77 = 1. The weaker 
characteristic scattering process is denoted by the dotted ar- 
row. Spins are indicated as block arrows. 



In the actual case, there will always be a region near 
normal incidence where the refiection amplitude will be- 
come linear (because it is antisymmetric). This means 
that 'very' far away ~ the inverse of the /c^-span of the 
region where r is linear, the previous result (|6| for the 
generic barrier should hold. 



C. Bi2Te3 (Avith hexagonal Avarping)^ 

If we are far away from the Dirac point, the surface 
band of Bi2Te3 exhibits hexagonal warping^ ^. The fol- 
lowing results hold when the STM bias maintains our 
observation energy in that regime. 



i. Step edge 1. FM direction 

The scattering wave-vector varies slowest around the 
equator (as in the 2DEG case), where Ao=diameter of 



2. Step edge 1. FK direction 

The wave- vectors vary slowest around the equator (as 
in the 2DEG case) and from the considerations of the 
cicular STI surface band above, we can conclude that 
there should be characteristic oscillations at 2krK decay- 
ing as x~^/^ (or 1/x for a 'perfect' refiector). In this 
case, because of the larger extent of the characteristic 
scattering region around the diameter, these oscillations 
may be strong and observable. Because of the presence 
of the linear band shape with larger spectral presence 
and refiection strengths near the corners, we can also 
observe LDOS oscillations from the corner^corner scat- 
tering processes, decaying as 1/x. Of course, from our 
simple calculation we cannot reliably predict which of 
the two processes discussed above have the stronger sig- 
nature. 




FIG. 5: (Color online) Scattering processes from a step edge 
on Bi2Te3 oriented perpendicular to the FK direction, in the 
energy range where the band exhibits hexagonal warping. 
Spins are indicated as block arrows. 



IV. ISOLATING CONTRIBUTIONS USING THE 
1-D FOURIER TRANSFORM 

The Fourier Transform of the LDOS data may be used 
to observe signatures from more than one set of scattering 
processes. For oscillations decaying as sm{Kx + 0)/x"^, 
scaling analysis tells us that the Fourier transform looks 
like F{k) ^ \k ^ K\^~'^, when k ~ ±K. Thus, one way 
to look for contributions to these oscillations would be 
to scan the 1-D Fourier transform of the LDOS (taken 



over the long-distance behavior) for features at the 'char- 
acteristic' scattering vectors. The FT near those points 
can then be fitted to the abovementioned power laws (or 
a logarithm, for the case of a 1/x decay) to recover the 
spatial decay power laws. 



V. CONCLUSION 

We have outlined a method to calculate the possible 
oscillatory power laws governing the decay of LDOS per- 
turbations next to a step edge or some such linear barrier 
on a surface, in the energy range when only electronic 
surface states are relevant and the surface band struc- 
ture is qualitatively known. To find these laws we need 
to identify the characteristic scattering regions (Figure 
[y, compute the scaling powers of the relevant quantities 

3) and from that obtain the possible oscillatory powers 

4). 
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